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ABSTRACT: The ability of unicellular green algal species such 
as Chlamydomonas reinhardtii to produce hydrogen gas via iron- 
hydrogenase is well known. However, the oxygen-sensitive 
hydrogenase is closely linked to the photosynthetic chain in 
such a way that hydrogen and oxygen production need to be 
separated temporally for sustained photo-production. Under 
illumination, sulfur-deprivation has been shown to accom- 
modate the production of hydrogen gas by partially- 
deactivating 0 2 evolution activity, leading to anaerobiosis 
in a sealed culture. As these facets are coupled, and the system 
complex, mathematical approaches potentially are of signifi- 
cant value since they may reveal improved or even optimal 
schemes for maximizing hydrogen production. Here, a 
mechanistic model of the system is constructed from 
consideration of the essential pathways and processes. The 
role of sulfur in photosynthesis (via PSII) and the storage and 
catabolism of endogenous substrate, and thus growth and 
decay of culture density, are explicitly modeled in order to 
describe and explore the complex interactions that lead to H 2 
production during sulfur-deprivation. As far as possible, 
functional forms and parameter values are determined or 
estimated from experimental data. The model is compared 
with published experimental studies and, encouragingly, 
qualitative agreement for trends in hydrogen yield and 
initiation time are found. It is then employed to probe 
optimal external sulfur and illumination conditions for 
hydrogen production, which are found to differ depending on 
whether a maximum yield of gas or initial production rate is 
required. The model constitutes a powerful theoretical tool 
for investigating novel sulfur cycling regimes that may 
ultimately be used to improve the commercial viability of 
hydrogen gas production from microorganisms. 
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Introduction 

Although the ability of the unicellular microorganism 
Chlamydomonas reinhardtii to photosynthetically produce 
hydrogen gas from water under illumination has been known 
for over 60 years (Gaffron and Rubin, 1942), until recently it 
remained largely a biological curiosity as hydrogen producing 
iron-hydrogenase is inhibited by oxygen co-produced from 
the photosynthetic pathway under normal illumination and 
nutrient conditions (Benemann et al, 1973; Ghirardi et al., 
1997, 2000). Thus photosynthetic growth and hydrogen 
production are incompatible and need to be spatially or 
temporally separated in order to achieve significant hydrogen 
production. Melis et al. (2000) proposed a groundbreaking 
two-stage process for temporally separating the hydrogen and 
oxygen components of the photosynthetic pathway: cells are 
grown as normal in a sulfur-replete media and then in a 
second non-growth stage, partial deactivation of the oxygen- 
evolving photosystem II (PSII) occurs in response to sulfur- 
deprivation. In essence, during water splitting in PSII, the 
sulfur- rich reaction-center Dl proteins are damaged and 
need to be replaced (Mattoo and Edelman, 1987). In the 
absence of sulfur, Dl protein biosynthesis is impeded and the 
PSII repair cycle is blocked (Wykoff et al., 1998), leading to a 
reduction in oxygen production to a low level (Melis 
et al., 2000). Aerobic respiration and the light-dependent 
activity of photosystem I (PSI) are not directly affected by 
sulfur-deprivation (Cao et al., 2001; Davies et al., 1994; Melis 
et al, 2000; Zhang and Melis, 2002). After approximately 24 h 
under illumination, the rate of oxygen produced from 
photosynthesis is less than the rate of oxygen consumed by 
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respiration; in a sealed container, the cells consume dissolved 
oxygen in the medium and the culture becomes anaerobic 
(Ghirardi et al., 2000; Kosourov et al., 2002; Melis et al., 2000; 
Zhang et al, 2002). In addition, during this time electrons 
result from the catabolism of endogenous substrates such as 
protein and starch (e.g., Chochois et al, 2009; Fouchard 
et al, 2005; Posewitz et al., 2004), both of which have been 
shown to increase significantly in the initial stages of sulfur- 
deprivation before hydrogen is produced (Fouchard 
et al, 2005; Kosourov et al., 2002; Melis et al., 2000; Posewitz 
et al., 2004). These events cause morphological changes in the 
cells during hydrogen production (Zhang et al., 2002). 
During dark fermentation ethanol acts as an electron sink for 
any reducing equivalents produced, but ethanol is harmful to 
the cell (Kennedy et al, 1992). In the light, under sulfur- 
deprivation, the partially active respiratory chain does not 
suffice as an electron sink and nor does the Calvin cycle since 
Rubisco, a necessary sulfur- rich enzyme in carbon fixation, is 
broken down and not synthesized (White and Melis, 2006; 
Zhang et al., 2002). The oxygen sensitive iron-hydrogenase 
enzyme on the thylakoid membrane is activated under these 
conditions and steps in as a major electron sink, re-oxidizing 
potentially harmful electrons produced from both the PSIT 
dependent (via water splitting) and the PSII-independent 
(fermentation) pathways, yielding H 2 gas for around 100 h in 
the light (Fouchard et al., 2005; Happe et al, 2002; 
Hemschemeier et al, 2008; Kosourov et al., 2002; Melis 
et al., 2000). The catabolic PSII-independent pathway is 
thought to contribute 20% of the hydrogen production and 
the PSII-dependent pathway contributes 80% (Fouchard 
et al, 2005; Volgusheva et al., 2013). Substantial hydrogen 
production ceases after around 120-140 h of sulfur-depriva- 
tion, thought to be due to depletion of the endogenous 
substrate available for catabolism (see Melis, 2002). Hence, 
there is a metabolic transition between an aerobic state with 
photosynthetic growth and an anaerobic state characterized 
by fermentation, H 2 production and biomass reduction 
(Hemschemeier et al., 2008, see also Fig. 1). If sulfur is added 
to the culture once hydrogen production has ceased, the cells, 
and particularly PSII, can repair; cycles of oxygen production 
under S-sufficiency and H 2 production under S-deprivation 
can result (e.g., Ghirardi et al, 2000). 

The above description of the interplay between cellular 
processes is a simplification of very complex dynamics that 
whilst gaining general acceptance in the research community 
is subject to improvement (for recent reviews see Antal 
et al, 2010; Ghysels and Franck, 2010). Although the 
promising sulfur- deprivation protocol allows for significant 
hydrogen production, the efficiency of the two-stage process 
and the yields of hydrogen need to be improved to allow for 
commercial exploitation (see for example Das and Veziroglu, 
2008; Melis, 2002). Scoma et al. (2012) demonstrated 
hydrogen production from green algae from solar light for 
the first time, but found that light conditions and mixing had 
a large effect on the H 2 yield (see also Giannelli et al., 2009) 
(which is expected since the collective swimming behavior of 
such species is sensitive to light conditions, which in turn 



affects photosynthetic efficiency; Bees and Croze, 2010; 
Williams and Bees, 2011b). Furthermore, a large downtime 
arises due to sulfur-cycling between anaerobic sulfur- 
deprived hydrogen production and aerobic, sulfur- replete 
recovery periods. In order to advance beyond the standard 
two-stage process it is first necessary to understand the 
system within the limits of this procedure. 

Strategies for the optimization of hydrogen gas production 
via the two-stage process can be designed and tested using 
dynamical models to represent the main pathways and 
processes of the system. To this end, we construct a simple 
mechanical, mathematical model of an algal culture that can 
describe sulfur-deprived hydrogen production in C, reinhardtii 
from a careful consideration of the biology and biochemistry, 
including important feedback pathways. The model is general 
in the sense that it captures both sulfur-deprived and sulfur- 
replete conditions. Beyond non-mechanistic approaches (Jo 
et al., 2006; Jorquera et al., 2008), there are two mechanistic 
models of aspects of the algal system under these conditions. 
Park and Moon (2007) constructed three separate state models 
of the biochemical photosynthetic processes involved in 
hydrogen production and specifically modeled eight primary 
metabolites. The release of hydrogen gas and the effects of 
illumination were explicitly modeled, but the role of 
endogenous substrates was omitted. Furthermore, the model 
is a discrete, multi-state model rather than a continuous 
formulation, and parameters values were difficult to identify. 
Fouchard et al. (2009) improved upon this approach by 
formulating a continuous description of the role of sulfur and 
light limitation in photosynthetic growth and anaerobiosis 
under general conditions and applied their model to the case of 
sulfur-deprivation, but the model stopped short of modeling 
the production of hydrogen gas. Model validation and 
optimization were considered by Degrenne et al. (201 1). 

In this study we shall improve upon previous work by 
modeling the principal mechanisms for the whole hydrogen 
production system, including feedback between sulfur uptake, 
photosynthetic growth, endogenous substrate, and the release of 
H 2 gas. There are elements that are modeled in a similar fashion 
to Fouchard et al. (2009). In particular, both intra- and extra- 
cellular sulfur are considered and we describe the uptake of 
external sulfur using a modified Monod formulation (Monod, 
1949) and illumination and photosynthetic activity are 
dynamically coupled, since it is well known that culture 
growth has an effect on the light available for photosynthesis. We 
describe the effects of sulfur- deprivation on the rate of 
photosynthesis using a similar modified- Droop relationship 
(Droop, 1968, 1979), and the use of sulfur in PSII repair, and the 
release of oxygen from PSII and its consumption in respiration 
are also included. But, significantly, the current approach 
extends and refines previous work in a number of ways. Firsdy, 
we model explicitly and mechanistically the initial storage and 
subsequent catabolism of endogenous substrate: protein 
breakdown in particular is important due to the release of 
small amounts of sulfur that can permit residual PSII activity, a 
key source of electrons for H 2 production (e.g., Fouchard 
et al., 2005; Melis et al., 2000). We model substrate storage as 
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Figure 1 . Schematics of the intracellular processes and pathways that occur under normal, sulfur-replete conditions (panel a) and during sulfur-deprivation (panel b). Light 
gray arrows and text indicate an inactive pathway/process. In panel (a), sufficient sulfur levels allow maximal PSII repair. Electron flow (dashed arrows) from PSII to PSI leads to ATP 
synthetase and oxygen production that inhibits the activity of the iron-hydrogenase (thick black line), where the Calvin cycle is active. Under sulfur-deprivation (panel b), PSII activity 
decreases, fermentation begins (releasing minimal quantities of sulfur and electrons) and low Calvin cycle activity, caused by Rubisco depletion, activates the iron-hydrogenase 
under anaerobic conditions. 
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dependent on the illuminated and S-dependent photosynthetic 
pathway (since proteins and starch are made via the Calvin 
cycle) and substrate breakdown (fermentation) as an emergency 
response to anaerobiosis, which also provides electrons to the 
hydrogenase. Culture growth can then be modeled as a function 
of endogenous substrate. These aspects differ from Degrenne 
et al. (2011) and Fouchard et al. (2009), since in these articles, 
changes in biomass are partitioned into growth and starch 
accumulation and these processes are not modeled indepen- 
dendy (both depend in the same way on photosynthetic rate), 
protein dynamics and fermentation are not modeled explicidy, 
and feedback between substrate catabolism and sulfur release 
for PSII is not incorporated. The model could not adequately 
capture observed starch accumulation dynamics (Degrenne 
et al., 2011). The new description also deviates from the 
previous culture growth models in that it allows for both culture 
growth and biomass reduction under nutrient limitation (as 
shown experimentally in Zhang et al., 2002). Furthermore, it 
provides feedback pathways between growth, substrate catabo- 
lism and S-dependent photosynthesis. And finally, hydrogen 
production is modeled explicitly as a system output that is 
dependent on light and electron donation via both the PSII- 
independent and PSII-dependent pathways and is inhibited by 
oxygen within the culture (Ghirardi et al., 1997, 2000). 

The model presented here consists of a set of coupled 
ordinary differential equations driven by evolving culture 
conditions (see Williams, 2009). In the following sections, the 
model is constructed from a mechanistic perspective and the 
solutions explored numerically. As in previous publications 
in this area, parameter estimation and the determination of 
functional forms were considerable challenges. Our objective 
was to produce a robust mechanistic model that exhibits 
the same qualitative trends as observed in experiments, 
rather than to refine parameter values arbitrarily to obtain 
quantitative agreement. Parameter values or estimated 
ranges were obtained from published experimental studies 
(see "Supplementary Material") and the model was then 
employed to probe the system subject to the constraints of 
the two-stage process outlined above. Model results are 
compared with published experimental data, and optimal 
external sulfur and illumination conditions are determined. 
In a subsequent paper, novel sulfur- cycling strategies will be 
explored for optimizing hydrogen production outside the 
confines of the two-stage process. 

Model Assumptions and Descriptions 
Model Formulation 

The mechanistic model that we shall develop consists of a set 
of mass balance equations that represent three stages of 
hydrogen production: normal photosynthetic growth 
(Fig. la), activity under sulfur- deprivation and subsequent 
hydrogen production (Fig. lb). We model an asynchronous 
cell population in a sealed, cubical container, purged of 
oxygen at t = 0 and filled with 1 L of culture, with illumination 
of 300|jiEm~ 2 at two sides. Cell division and changes in 



individual cell size are combined into one variable, the cell 
volume fraction 0 < A < 1 . We assume that oxygen diffusion 
across the cell wall is rapid, and thus internal and external 
oxygen concentrations can be described by one variable a>, in 
pM. External and internal sulfur, S and s, respectively, are not 
combined as they have distinct dynamics with active sulfur 
transport across the cell wall (Yildiz et al, 1994). Endogenous 
substrate, e, and protein p are also modeled explicidy. The 
variables s, e, and p are intracellular concentrations (or quota, 
Droop, 1968, 1979), in pM: concentrations within the 
suspension in the bioreactor are given by, for example, pA . S is 
modeled as concentration in the suspension medium (pM; 
the concentration of external sulfur in the suspension is thus 
S(l — A)). Hydrogen gas concentration, h, is modeled as a 
product in units of mL/L of culture, where we assume the 
same experimental and altitude conditions as Kosourov et al. 
(2002) such that 1 mL H 2 = 33 pmolH 2 . 

The Effects of Culture Density on Light Availability 

Cell volume fraction affects the optical density of the culture 
and thus the rate of photosynthesis. We model light intensity 
using the Beer-Lambert law (see Duysens, 1956), assuming 
that the cells are homogeneous and transmit light equally in 
all directions. The small effects of multiple scattering are 
neglected. Illumination is from the side, at x = 0, where 
0 < x < b w measures the distance from the container edge to 
the light source of intensity / 0 (b w is bioreactor width). 
Assuming uniform cell concentration, n(x) = n 0 , light 
intensity is given by I(x) = 7 0 exp(A; c hi«o-" : )» where for 
simplicity, the absorbance of the medium is assumed 
negligible and A; cm is the absorbance of the cells. Furthermore, 
we assume that the culture is well-mixed so that averaging 
over the width is the same as averaging over time: 



(I) 



1 f° 

, = (I) x = i- h exp(k chl n 0 x) dx 

K J-b„ 



(1) 



where (/) ( and (I) x indicate time and space averages of /, 
respectively. These assumptions are incorrect if mixing is 
weak or swimming induced bioconvection results (Bees and 
Croze, 2010; Williams and Bees, 2011b). Furthermore, there 
is a light intensity, / sat at which the photosynthetic rate 
saturates (Leverenz et al., 1990). Hence, imposing a Heaviside 
function and integrating (Supplementary Material; Williams, 
2009), the dimensionless usable light, L(A), is 



L(A) = 



D C A 



exp 



-In 



^sat In 



x H\ 1 



D C A 



In 



-exp(-Dc^) 



D C A 



(2) 



where we define D c = fr w ^chi/^ceii- Here, L(A) has been 
normalized with I el and the light intensities are non- 
dimensionalized with L el (e.g., / 0 = h/L e i), a standard value 
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Table II. Expected range of non-dimensional parameter values, calculated from Table I. 



Dimensionless parameter 


Definition 


Description 


Standard value 


Range 






Scaled ratio of sulfur uptake rates 


1 £LA \j 1 ffi 

1.64 X 1U 


o cia i n4 a c », i n' 
H.y4 X 1U —5.5 X IU 


A 2 




Scaled ratio of sulfur uptake rates 


159 0 


Q (SASH 
y. 3y— t)4oU 


B 




Ratio of S uptake rates 


6.65 


a an i a ^ 


D c 


v„i 


Scaled measure of absorption 


26,300 


z X 1U — o X 1U 




El 


Fraction of electrons from PSII-dependent path 


0.75 


0.7—0.8 




In, 
i,i 


Normalized light saturation 


24 8 


20—30 0 


J o 


fo 

1,1 


Normalized light intensity at source 


99 2 




ir 
A 2 


h 
h 


Scaled protein breakdown rate 


1 95 




ir 
K 3 


Ml 


Scaled measure of rate of p production 




U.UUD /—ZD. 4 


^5 


h 

h 


Scaled respiration rate 


i II „ 1 n<* 

o.44 X 1U 


D.4fl X 1U —/A/ X 1U 


^6 


h 
hx 


Scaled photosynthesis rate 


1.2 x 10 


7 an v 1 1 Q7 v in' 

/.jy x iu — i.y/ x iu 


^G 


h 


Scaled protein gradient 


3.52 


u.oyy— 3.y 


P H 




Non-dimensional reciprocal of pj, 


0.560 


n aon "7 "7Q 

u. dou— z.zy 


P K 


h 


Scaled protein required for survival p r 


0.292 


U.U43/— U.DDD 


Kg 


>-.-y.y 

h 


Scaled growth rate 


1 56 


0 82—2 39 


#D 


'■decay 
*1 


Scaled decay rate 


0.129 


n nil i n 


Si 


£l_ 

«n 


Ratio of s required for Calvin cycle to normal s 


0.5 


0.2—1.0 


Sh 


Sa 
«h 


Non-dimensional reciprocal of $h 


6.0 


4.0—12.0 


V L 


vo 2 


Scaled oxygen mass transfer 


9.12 


0.665-13.3 


Yo 




Scaled protein switch p 0 


1.94 


0.542-2.58 


Yi 


*3 


Scaled protein switch p l 


1.91 


0.516-3.18 


Yi 


*3 


Scaled protein switch p 2 


2.23 


0.647-3.11 


[2 1 


0», 

X 


Scaled oxygen switch coj 


0.0047 


0.0025-0.01 


S2 2 


m 

X 


Scaled oxygen switch co 2 


0.103 


0.0433-0.195 



S, hydrogen, h, rate of hydrogen production, d/i/df, and 
time. 

Inside the cell, sulfur is used in replacing photo-damaged 
PSII (termed "repair" herein) and in making other proteins. 
We assume the use of sulfur for PSII repair is linearly 
dependent on light, due to photo-damage, and available 
sulfur. A Heaviside switch function H PSU denotes that above a 
critical concentration of internal sulfur, s n , photosynthetic 
activity is not affected by s concentration, and photosynthesis 
and thus PSII repair occurs at a constant rate: 

repair(s, A,L(A)) = -kiA(sH m i{s n s) + s 0 H V sii(s-s 0 ))L(A) 

(5) 

This relationship is analogous to the modified Droop 
formulation with a switch function employed by Degrenne 
et al. (201 1) and Fouchard et al. (2009) in which PSII activity 
drops off rapidly once sulfur falls below the critical quota 
value. Encouragingly, the corresponding curve of photosyn- 
thetic activity as a function of s agrees with the experimental 
measurements of Fouchard et al. (2009) (not shown). 

Intracellular protein concentration p is a large component 
of endogenous substrate and can act as a sulfur store: during 
anaerobic fermentation protein is catabolized to release sulfur 
(Melis et al, 2000). We model this sulfur source as dependent 
on available protein and oxygen levels using a switch function 
to specify that fermentation only occurs during anaerobiosis, 
co < &> p (Happe et al, 2002). A non-consumable base level of 
protein, p„ necessary for cell survival is also modeled (shown 
experimentally in Kosourov et al, 2002) 

protein breakdown(p, &>, A) = k 2 (p~p t )AH Ferment (co v ~oo) 

(6) 



employed in Kosourov et al. (2002) subject to which other 
parameters are measured and inferred (tildes have been 
dropped, see Table I and Supplementary Material). 

Sulfur Kinetics 

We employ data from Yildiz et al. (1994), to model the uptake of 
external sulfur into the cells from the media: sulfur uptake is 
dependent on both external and internal sulfur concentrations 
(uptake rate varied between sulfur-starved and normal cells; also 
shown experimentally in Fouchard et al., 2009). This leads to a 
modified Monod formulation for the total sulfur uptake for cell 
volume fraction A in which the Michaelis-Menten uptake rate 
under normal, s-replete conditions, ce(s), and the half saturation 
value, P(s), are in this case sulfur dependent functions: 

uptake (S,s,A) = (3) 

Assuming that 5 = 0 in the starved cells and that s = s n , the 
"normal" amount of sulfur, within an unstarved cell then we fit 
a(s) = aexp(-Gs/s n ), fi(s) = b\ + b 2 (s/s n ) to the data in 
Yildiz et al. (1994), where a is the maximum uptake rate of 
external sulfur (values shown in Table I). Thus for total external 
sulfur in the media, S(l — A), we obtain 

aA ^> +F ( s , h M t] (4) 



v S Input 

uptake 



where F is an arbitrary addition of external sulfur to 
the bioreactor, which may depend on external sulfur, 
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Figure 2. A diagram of the transport of electrons from the PSII-dependent pathway and the PSII-independent pathway to PSI and on to either the iron-hydrogenase, 
fermentation products or the Calvin cycle. Where the electrons end up is decided by the oxygen and sulfur dependence, as indicated. 



Protein is produced under normal conditions, combining 
sulfur with carbon skeletons produced from the photosyn- 
thetically dependent Calvin cycle. Thus we model protein 
production linearly on sulfur availability and light intensity 
up to a certain concentration of s, s n , using ffpsn normalized 
with s n to stipulate that given sufficient sulfur, photosynthetic 
activity is constant (see Equation 5). We assume that the 
cell can use one of the Calvin cycle, ethanol production or 
H 2 production as an electron sink at any one time (see Fig. 2). 
This assumption is realized by using a switch function 
-Heal™ (s—Si) to stipulate that protein is only produced 
when sufficient sulfur (s>s 1 ), thus sufficient Rubisco, 
allows the Calvin cycle to function (White and Melis, 
2006; Zhang et al, 2002) (see Disussions and Conclusions 
section). Thus the model of internal 5 andp concentrations is 



where /S indicates that 1 mol of protein contains /3 moles of 
sulfur. The PSII repair term does not appear elsewhere in the 
model, implying that sulfur used in PSII repair is not recycled. 
Note the growth terms in Equations (7) and (8), which arise as 
increasing cell volume alone reduces concentration. 

Oxygen Kinetics 

Under normal conditions PSII produces oxygen and 
respiration consumes oxygen. The relationship between 
PSII activity and sulfur is given in Equation (5) (for the rate of 
sulfur consumption) and is used here but with oxygenic 
photosynthetic rate constant k 6 (ixM0 2 h~ 1 ) and the non- 
dimensional PSII- switch. Respiration rate remains relatively 
unaffected by sulfur-deprivation (for t<70h; Melis 
et al., 2000) and thus is modeled as constant when oxygen 
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- — L(A)(sH VS ll(S n ~S)+ 5 n Hpsn(s - 5 n ))tf C alvin(-S -Si)- — — 

protein production growth/decay 



d ^3 

f t = ~k 2 (p ~ pr ) ^Ferment Op - O)) + — L(A) (s/fpsn (s n ~ s) + 5 n ^ PS Il(s - 5 n ) )H Ca lvm (s - Si) 
* ' $n 



protein breakdown v 

protein production , . 

p dA ( ' 

growth / decay 



Biotechnology and Bioengineering, Vol. Ill, No. 2, 2013 



is sufficient, a>> u> u but decreases linearly when oxygen is 
sparse, aXco^ Although the bio-reactor is sealed, we 
stipulate that 0 2 can leave the system when the culture is 
oxygen saturated and cannot reenter. Thus 



Hydrogen Production 

Hydrogen production is modeled as dependent on the scaled 
sum of electrons coming from endogenous substrate 
catabolism (the PSII-independent pathway, dependent on 



f = Ak b (j-H PSII (s n -s) + ffp Sn (5 - s n )^jL(A) 



photosynthesis (9) 

- Ak 5 (<MH Resp (a>i - co) + &>itf Resp ((W - a>i)) - v 0l {m - x)Hl oss (oj - x) 

V ^ / v ' 

respiration super— saturation loss 



Growth and Decay of Cell Culture 

Changes in cell volume fraction can occur due to both 
photosynthetic biomass production and substrate break- 
down. Since endogenous substrate is explicitly modeled, we 
can incorporate changes in culture concentration dependent 
on endogenous substrate rather than on the rate of 
photosynthesis (as modeled in Fouchard et al., 2009) to 
allow for both growth and decay of the culture (Zhang 
et al., 2002). Under sulfur sufficiency the culture density 
increases as endogenous substrate increases, whereas under 



dh 
dt 



k 4 A H Sensi , ivity (tti2 — <w) 



0 2 sensitivity 



p) and from the residual level of the PSII activity (the PSII- 
dependent pathway, dependent on 5) (Fouchard et al, 2005; 
Happe et al, 2002). These pathways are assumed indepen- 
dent and, as for growth, we assume protein and starch 
catabolism act in the same way and base our model on 
protein concentration alone. The Calvin cycle switch 
function ifcaivin( s i — s ) is used to stipulate that the Calvin 
cycle also needs to be inactive for the hydrogenase to function 
as the electron sink (see Fig. 2), and an oxygen sensitive 
switch is employed to reflect the dependence of iron 
hydrogenase on anaerobic conditions. Thus 



£ L - +(1-£ L ) ( ^ ) 
v _j^ ^ ph 

\PSII-dep PSII-indep / 



L(A) H Calvin (s, - s) 

" v ' 

e" pathway 



(11) 



nutrient deprivation substrate breaks down and the culture 
density decreases. Changes in endogenous substrate are 
largely due to protein and starch, both of which initially 
increase and then decrease during H 2 production (e.g., 
Degrenne et al., 2011; Fouchard et al., 2005; Kosourov 
et al., 2002). We assume that protein and starch are 
sufficiently correlated to allow modeling them as one entity, 
and thus we model growth rate as explicitly dependent on 
protein (e — p from here on). 

The growth function is chosen so that the growth and 
decay rates are constant above and below, respectively, critical 
levels of protein (p 2 and pi, respectively, see Table I for 
details); and there is some linear transition between the two 
states, giving an "s" shaped function (using a smoothed 
version had no qualitative effect). Light dependence is 
modeled explicitly in protein production and thus is implicit 
in the growth term. Thus 

cL4 
dT 



A 



r ap H G2 (p - p 2 ) + MP P " ) H G2 (p 2 - p)H G1 

(P2 ~ p0) 



(P-Pl)~ fdecay^GlOl ~ P) 



(10) 



where r exp and r decay are the maximum growth and decay 
rates, respectively. 



where E L is the fraction of electrons from the PSII-dependent 
pathway under total sulfur- deprivation (see Table I and 
Supplementary material for details). 

Thus mass balance equations (4) and (7)— (1 1) make up 
the standard model. Parameters definitions, values, and 
references are summarized in Table I. Parameters are taken 
from the literature where possible or else estimated using 
available relevant data (see Supplementary Material for full 
details). The model is non-dimensionalized using t = k\t, 
S = S/b 2 , 1 = s/s n , p = k 2 p/ks, co = ao/ x> and h = k\hjk±. 
The scaling for time is chosen so that one non-dimensional 
time unit corresponds to approximately 1 day. The sulfur 
scaling is chosen so that 5=1 initially (under normal 
conditions). The non-dimensional standard model equa- 
tions and parameters are shown in Appendix A (tildes are 
dropped from here on). 

Results 

To illustrate the model dynamics the set of differential 
equations were solved numerically using Matlab 7.0 software 
(R2007) with the robust implicit scheme "ode 15s" (employ- 
ing a modified backward Euler method; Shampine and 
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Reichelt, 1997). The numerical method was verified with a 
known solution for a simplified and linearized version of the 
model. The initial conditions were chosen to be representa- 
tive of experimental conditions in Kosourov et al. (2002): 
5=1, a>= 0, p = 2.23, A = 2.25 x 10" 3 , h = 0, and initial 
external sulfur concentration So varied between model runs. 
For numerical simulations, approximations to continuous 
functions are preferred over discontinuous (Heaviside) 
switches (as in Degrenne et al, 2011). Thus hyperbolic 
tangent switches were thus employed: 

H switch (F c -F) = i (1 + tanh(g(P c -F))) (12) 

This function of F varies rapidly from 1 to 0 around the 
critical value F c for a large value of the parameter g; it was 
increased to a value beyond which it did not significantly 
affect model output (Williams, 2009). 

Results with a large cell concentration in a sulfur- replete 
medium in sealed conditions are shown in Figure 3. There is a 
rapid increase in cell volume fraction in the first 3 days, with a 
doubling time of 22-23 h, compared to 9.4-18.6 h calculated 
from experimental data (Fischer et al., 2006; Jo et al, 2006). 
After 3 days, light limitation decreases oxygen production 
from PSII, resulting in anaerobic fermentation. Hydrogen 
production is not observed as sufficient sulfur is available for 
the Calvin cycle (via Rubisco) to act as the electron sink. 
Fermentation causes cell volume fraction to decrease, 
decreasing light limitation and, since sulfur is available, 



PSII activity increases and the system subsequently becomes 
aerobic and a period of protein production and growth 
follow. Thereafter, oscillations in s, p, a>, and A are found 
(period around 97 h), with no hydrogen produced. These 
results are consistent with Zhang et al. (2002), in which a 
concentrated culture in a sealed container became anaerobic 
as cell density increased, but only inactive hydrogenase was 
found. A bioreactor culture could be continuously diluted to 
optimize growth and avoid over-densification (see also 
Fouchard et al., 2009). 

Figure 4 shows the model results for a culture suspended in 
a sulfur- free media, So = 0 (xM, at t= 0. The model is run for 
approximately 10 days. Internal sulfur immediately starts to 
decrease while cell volume fraction increases initially as sulfur 
is still sufficient for Calvin cycle activity and growth, s > S t . 
When s falls below S : growth slows down as the Calvin cycle 
becomes inactive due to a lack of sulfur and thus protein is 
not produced. As s decreases further, the oxygenic photosyn- 
thetic rate falls below the respiration rate and a period of 
anaerobiosis begins after approximately 1 day. Since the 
Calvin cycle is also inactive under s-deprivation, hydrogen 
production now commences, and fermentative protein 
breakdown begins, resulting in release of small amounts of 
internal sulfur, p and A decrease during this H 2 production 
phase due to catabolism of endogenous substrate, p reaches 
P R , the base level of protein needed for cell survival, between 2 
and 4 days, but protein breakdown continues to supply 
electrons and sulfur to the photosynthetic pathway because 
the shrinking cell volume fraction causes oscillations in p 
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around P R (as A decreases, cellular protein concentration p 
increases transiently; total protein in the culture, pA, 
monotonically decreases). The initial hydrogen production 
rate is rapid but decreases significantly at around 6 days, when 
internal sulfur has run out, PSII activity stops, endogenous 
substrate is low and only minimal amounts of hydrogen are 
now produced from the PSII-independent pathway After 2 
more days hydrogen production stops and the cells continue 
to shrink. The final yield of gas after ten days is 106 mL H 2 /L 
culture, and after 140 h (f=5.74) is 103mLH 2 /L culture 
compared to 7 1 . 7 mL H 2 /L culture in 1 40 h in Kosourov et al. 
(2002). Results for s, a> and A are qualitatively similar to 
model results by Fouchard et al. (2009). 

Optimizing H 2 Yield: Varying Initial External Sulfur, S 0 

Re-suspending the cells in media with minimal rather than 
zero concentrations of external sulfur has been shown to 
increase the total yield of hydrogen gas (Kosourov et al, 2002; 
Zhang et al, 2002). Figure 5 shows model results for different 
initial concentrations of external sulfur, S 0 . For S 0 > 0 internal 
sulfur and protein decrease slower than when S 0 > 0, leading 
to higher culture density. Increased oxygen combined with a 
later decay in p and s leads to a later onset of anaerobiosis and 
a delay between this onset and hydrogen production. For 
S 0 =3.45 (50 u-M), yields of hydrogen gas are significantly 
larger than for S 0 >0 at f = 10 (/i = 237mLH 2 /L vs. h= 106 
mL H 2 /L, respectively), and production begins later (f = 45.2 
h when So = 3.45 and f = 36.4 h when S 0 = 1.725). Figure 6 



shows these results in detail: increasing S 0 from zero to 
S 0 = 6.9 delays the onset of H 2 production and increases yield 
at t = 10 but, as S 0 is increased further, yields decrease until 
hydrogen is not produced in this time frame. The optimal S 0 
for H 2 yield at t = 10 is S 0 = 6.19 (89.8 uM) with 6 = 246 
mL H 2 /L culture. 

The average initial rate of H 2 production over the first 15 h 
of production is calculated per unit of cell volume fraction 
using 

H 2 rate = ^ w (13) 

where T H is the onset time of hydrogen production, T l = 0.6 is 
the scaled initial time-period of production considered, and 
A is averaged over the initial hydrogen production period. 
Figure 7 shows a slight increase in the initial rate of hydrogen 
production per cell volume fraction as S 0 is increased from 
zero up to S 0 « 1.25, and thereafter the H 2 production rate 
decreases and reaches very low levels at S 0 = 6. Thus there is 
an optimal initial value for external S 0 for improving the rate 
of hydrogen produced per cell (S 0 ~ 1), which is different to 
the optimal for improving yield at time t = 10 (S 0 = 6.19). 

Optimizing H 2 Yield: Varying Light Intensity, l 0 

Varying the illumination conditions of the algal culture can 
have an effect on the yield of hydrogen gas (e.g., Degrenne 
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Figure 5. Results for the model with standard parameter values in Table I, with initial conditions of external sulfur of S 0 = 0 (solid lines), S 0 = 1.725 (dotted lines), and S 0 = 3.45 
(dashed lines). These correspond to 0, 25, and 50 p.M, respectively. 



et al., 2011; Kim et al, 2006). Figure 8 shows results for three 
values of the light intensity, 7 0 when S 0 = 3.45 (50 u,M). For 
7 0 = 49.6 (half of the standard value), slower growth and, 
hence, slower s usage delay onset of H 2 production. A smaller 
cell volume fraction combined with reduced activity of PSII- 
dependent activity result in a decreased H 2 yield, as expected. 
When I 0 is doubled from the standard value, I 0 = 198.4, rapid 
growth leads to higher cell volume fraction and faster sulfur 
usage compared to the standard case. However, the resulting 
increase in oxygen production causes the system to become 
anaerobic and hence produce H 2 at approximately the same 
time. Perhaps surprisingly, hydrogen production stops 



sooner and the yield is significantly reduced with higher I 0 
even though the cell density and light available for the PSII- 
dependent electron pathway have increased. This is due to 
increased PSII photo-damage causing a more rapid decline in 
internal sulfur, which limits PSII-dependent electron dona- 
tion. Thus, there is an optimal light intensity to maximize 
hydrogen yield within a given time (as shown experimentally 
by Kim et al. (2006)): for the model presented here, for 
fend= 10 the optimal light intensity is I 0 = 146.5 u,molm _2 
s" 1 if S 0 = 0 and I 0 = 340 |xmolm" 2 s _1 if 5 0 = 3.45 (50 uM: 
see Fig. 9). Thus the initial sulfur concentration has an effect 
on the optimal light intensity. Decreasing the cellular 
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Figure 6. Hydrogen and cell volume fraction curves for the model with standard parameter values in Table I and with initial conditions S o = 0 (solid lines), 3.45 (dashed lines), 6.9 
(dot-dashed lines), 13.8 (dotted lines), and 20.7 (thick dashed lines) in non-dimensional units. 
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Discussion and Conclusions 
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Figure 7. Initial rates of hydrogen production (in the first 15 h) plotted against the 
initial amount of external sulfur for the standard parameter values in Table I. 



absorption coefficient, D c , provides more light on average to 
each cell and thus also results in a greater hydrogen yield: for 
I 0 = 300 u-molm -2 s _1 , decreasing D c increases yield. How- 
ever, for large light intensity 7 0 and small D c , the effects of 
photo-damage cause an overall decrease in yield. 



A simple mechanistic model has been constructed to describe 
sulfur-deprived hydrogen production in green algae. By 
modeling mechanistically, we have significantly simplified 
this complex system to just six variables. Key features of the 
model, including sulfur-dependent photosynthesis, growth, 
changes in endogenous substrate, and hydrogen gas release, 
have been incorporated. Solutions were obtained for the 
standard values of the parameters and with a range of initial 
conditions. 

The experimental studies of Kosourov et al. (2002), Melis 
et al. (2000), and Zhang et al. (2002) guided the construction 
of the model and some parameters in the growth and 
hydrogen functions, in particular, were extrapolated from the 
experiments therein. For example, the hydrogen production 
rate constant k 4 was taken from measured data in Kosourov 
et al. (2002) for the case of S o = 0|jlM, thus the simulated 
hydrogen dynamics match the experimental data for So = 0 
|xM reasonably well, as expected. However, the model was 
not fit to the data, and the hydrogen dynamics for different 
initial external sulfur and illumination conditions can still be 
compared independently with experimental data in order to 
test whether the model correctly captures the system 
dynamics under different conditions. Better independent 
measurements of the parameters, rather than fitting, should 
be the focus of future research efforts. Encouragingly, good 
qualitative and quantitative agreement was obtained between 
experimental results and model simulations for H 2 yield for 
different initial external sulfur concentrations, So'- after 140 h 
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Figure 9. Hydrogen yield at f=10 as a function of dimensional / 0 when Sb = 0 
(solid line) and S 0 = 3.45 (dashed line) for the standard parameter values in Table I. 



if S 0 = 25 |jlM or S 0 — 50 |xM, the model predicts yields of 
/2=168mLH 2 /L culture and /2 = 213mLH 2 /L culture, 
respectively, in good agreement with Kosourov et al. 
(2002) (h = 127 and h= 159mLH 2 /L culture, respectively). 
The optimal S 0 for maximum hydrogen output over a fixed 
period was found to be a dynamic balance between high 
culture density, light limitation, and production start time. 
Hydrogen production onset time also corresponded approx- 
imately to experimental results: for So = 25 |xM and So = 50 
u,M, f=36.4, and f=45.2h, respectively, compared to 
f= 43-49 h in Kosourov et al. (2002). In simulation results, 
hydrogen production began almost as soon as the system 
became anaerobic when S = 0 (xM, as in Zhang et al. (2002), 
but Kosourov et al. (2002) found a slight delay between onset 
of anaerobiosis and hydrogen production. This delay was 
predicted by our model for S > 0 uM, due to slower sulfur 
decay causing an extended period of Calvin cycle activity. 

The initial rate of hydrogen production per cell was also 
investigated, and we found that it increased slightly then 
decreased substantially as S 0 increased. The relatively constant 
production rate per cell is consistent with experimental 
observations from Degrenne et al. (2011) and Zhang et al. 
(2002) (where rate is per gram of biomass), but inconsistent 
with Kosourov et al. (2002), who found an increase in initial 
H 2 production rate per mole of chlorophyll for S 0 = 25 jjlM 
compared to S 0 = 0. We attribute increased hydrogen yield for 
So « 50 u,M to increased cell volume fraction, as found 
experimentally by Zhang et al. (2002), rather than increased 
production rate per cell as proposed by Kosourov et al. 
(2002). The decrease in the initial rate for So > 43.5 jjlM 
found from our model is also consistent with the trends 
found by Kosourov et al. (2002) and Zhang et al. (2002) for 
So > 25 u,M. Likewise, we attribute corresponding decreases 



in H 2 yield to increased light limitation counteracting further 
increases in A when S 0 is large. The optimal sulfur 
concentration for maximizing this H 2 production rate 
(approximately 0 < S 0 < 29 (xM) was found to be different 
from the optimal sulfur for increasing overall yield (S 0 = 89.9 
uM). Thus methods of optimization of the hydrogen 
production system depend on whether maximum cell activity 
or maximum H 2 output per culture is required. 

Model simulations for changes in illumination are 
consistent with experimental data from Hahn et al. (2004) 
and Kim et al. (2006): increasing the light intensity / 0 can 
significantly increase yields up to an optimal value due to 
earlier onset of production and increased culture density and 
PSII-dependent electron flow. However, increasing 7 0 beyond 
the optimal value decreases H 2 yields due to increased photo- 
damage, as in Kim et al. (2006). Simulation results predict an 
optimal light intensity for total H 2 output of / 0 = 146.5 
u,molm _2 s _1 for So = 0|xM or I 0 = 340 jjimolm -2 s _1 for 
S 0 = 50 uM with illumination from both sides, which are of 
the same order as those predicted by Park and Moon (2007) 
(238(xEm" 2 s _1 ) and Kim et al. (2006) (200 |xEm" 2 s _1 ). 

Using the model of Degrenne et al. (2011), Fouchard et al. 
(2009) found qualitatively similar results: they predicted that 
a high hydrogen yield would require high external sulfur and 
light irradiance. Experimental data supported this conclu- 
sion. However, in those studies H 2 gas production was not 
explicitly modelled but was extrapolated from biomass and 
starch concentrations. We find that higher yields of H 2 are 
found for higher cell volume fraction: for S 0 = 0, h = 247 
mL H 2 /L with A 0 = 0.0045 and 106mLH 2 /L culture when 
A 0 = 0.00225, which supports the hypothesis that one may 
optimize H 2 yield by maximizing biomass. However, we 
caution that for sufficiently high initial sulfur and light 
conditions our model predicts diminished H 2 yields due to 
over-concentrated cultures and photo-damage. 

Melis (2002), Melis (2009), and Polle et al. (2002) 
suggested that truncating the chlorophyll antenna to decrease 
cellular absorbance (modelled as D c ) decreases wasted light 
and increases photosynthetic activity, which may increase the 
hydrogen yield. Model results also suggest that decreasing 
absorbance could optimize H 2 yield, provided that the light 
intensity 7 0 is not too high, or D c is not too low, otherwise 
yields decrease due to increased photo-damage (as for high 
light intensities in this model and in Kim et al. (2006) and 
Park and Moon (2007). 

To our knowledge, this is the first simple mechanistic 
model of sulfur-deprived hydrogen production to include 
feedback between sulfur, photosynthetic growth, endogenous 
substrate, and hydrogen production. Good qualitative 
agreement is found between model simulations and 
experimental results. In order to model such a complex 
system, key assumptions were made. The role of starch was 
not modelled independently; instead, endogenous substrate 
is representative of both protein and starch in order to 
capture the dynamical feedback between sulfur, photosyn- 
thetic growth and fermentation. This may be a reasonable 
approximation, but the two may be better modelled 
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separately, with growth a function of both. However, we do 
not expect this extension qualitatively to alter results. 

Additionally, a switch (H Ca ivin (s; Si)) was used to close the 
system and specify that H 2 -producing hydrogenase requires 
both anaerobiosis and an inactive Calvin cycle to function as 
an electron sink (e.g., Happe et al., 2002; Hemschemeier 
et al., 2008; White and Melis, 2006), so a sealed system with 
high culture density leads to anaerobiosis due to light 
limitations but no hydrogen is produced (in accordance with 
Zhang et al., 2002). In this study, the switch had little effect 
when initial external sulfur was minimal and it allowed the 
omission of the complex Calvin cycle and the interplay 
between electron sinks from the model. It may be revealing to 
explore further the explicit nature of the coupling between 
the hydrogenase and the Calvin cycle. 

To describe the suspension, the cultures were assumed to 
be perfectly mixed and cell swimming behaviour was not 
described. Biased swimming is known to induce hydrody- 
namic instabilities, resulting in non-uniform distributions of 
cells, called bioconvection, in tens of seconds on length scales 
of centimeters. This significantly affects light transmittance 
and thus photosynthesis (Bees and Croze, 2010; Williams and 
Bees, 2011a,b), and could have a substantial impact on H 2 
yield. All of these assumptions should be explored in future 
developments of the current model. 

As new data emerge, refinements of the parameter values 
and key mechanisms can be incorporated in the model. 
Perhaps more importantly, the current description is 
ideal for examining novel regimes for optimizing the total 
yield, or rates of production, of hydrogen gas produced 
under a range of sulfur-deprivation schemes. Such analysis 
may provide valuable insight into future commercialization 
of algal H 2 production and will be presented in a future 
article. 

Nomenclature 



a maximum rate of uptake of external sulfur 

(|xMh _1 ) 

Ai dimensionless ratio of sulfur uptake rates 

A 2 dimensionless ratio of sulfur uptake rates 

ATP adenosine triphosphate 

bi rate constants for sulfur uptake (i= 1, 

uM) 

fo w width of the bio-reactor (cm) 

B dimensionless ratio of S uptake rates 

D c dimensionless absorption measure 

e endogenous substrate (model variable) 

£ L fraction of electrons from PSII-dependent 

path 

F(S, h, dh/dt, t) input of external sulfur function (Equa- 
tion 4) 

g gradient of the smoothed Heaviside switch 

function 

G dimensionless scale factor in a(s) 

h hydrogen gas (model variable) (mLL -1 

culture) 



#"name" heaviside function to model process, 

"name" 

/ 0 dimensionless light intensity at the source 

/ sat dimensionless saturation level of light 

(I) x function for total light intensity averaged 

over width 

ki rate constant for PSII repair (h _1 ) 

k 2 rate constant for protein breakdown (h -1 ) 

k 3 rate constant for protein production (u,M 

h" 1 ) 

k 4 rate constant for hydrogen production 

(mLh -1 ) 

k 5 rate constant for oxygen consumption by 

respiration (h -1 ) 
k 6 rate constant for oxygen production from 

PSII (uMh" 1 ) 

fc chl measure of absorbance of the cells (cm 2 ) 

K 3 dimensionless measure of rate of protein 

production 

K 2 dimensionless protein breakdown rate 

K 5 dimensionless respiration rate 

K 6 dimensionless photosynthesis rate 

L(A) dimensionless usable light intensity 

(Equation 2) 

L el value of the usable light function in 

Kosourov et al. (2002) (|jimolm _2 s _1 ) 

n cell concentration (cells/mL) 

n 0 uniform cell concentration throughout the 

layer (cells/mL) 

NADP + nicotinamide adenine dinucleotide phos- 

phate 

NADPH reduced form of NADP+ 

p protein (model variable) (jjlM) 

p 0 protein level when growth is zero (u,M) 

pi protein below which maximum decay 

occurs (jjlM) 

p 2 protein required for maximum growth 

ph normalization of PSITindependent elec- 

tron pathway (u,M) 
p r basic protein needed for cell survival (|xM) 

P G dimensionless protein gradient 

P H dimensionless reciprocal of ph 

P R dimensionless protein required for survival 

PSI/PSII photosystem I/II 

r decay maximum rate for cell decay (h -1 ) 

r exp maximum growth rate (h _1 ) 

R u dimensionless decay rate 

R G dimensionless growth rate 

5 internal sulfur (model variable) (u,M) 

5 n normal level of sulfur in a cell (u,M) 

Si sulfur level above which Calvin cycle is 

active (|xM) 

5h normalization of PSII-dependent electron 

pathway (|xM) 
S external sulfur (model variable) (|xM) 

S 0 initial external sulfur (u,M) 
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Si 

S H 

t 

T 

Tu 
Ti 

v o 2 

^cell 
^cont 

v L 

a(s) 

P(s) 

P 

Yo 
Yi 
Yi 
A 

A T H 

X 

CO 
(Or 

co 2 

(Dp 

^2 1 

C2 2 



ratio of sulfur required for Calvin cycle 
compared to normal sulfur concentration 
dimensionless reciprocal of 5h 
time 

time at which total hydrogen yield h is 
output 

start time of hydrogen production 

scaled time period over which initial 

hydrogen production rate is measured 

oxygen mass transfer coefficient 

volume of a single cell (mL) 

volume of the container (mL) 

dimensionless oxygen mass transfer 

function for the maximum uptake rate of 

external sulfur (Equation 3) 

function for the half saturation value of S 

uptake (Equation 3) 

average moles of sulfur in one mole of 

protein 

dimensionless protein switch p 0 
dimensionless protein switch pi 
dimensionless protein switch p 2 
cell volume fraction (model variable) 
cell volume fraction when hydrogen pro- 
duction begins 

oxygen saturation in water (|xM) 
oxygen (model variable) (u,M) 
oxygen level required for full respiration 
(tiM) 

oxygen level required to inhibit H 2 pro- 
duction (|xM) 

oxygen level below which protein break- 
down occurs (|xM) 
dimensionless oxygen switch co x 
dimensionless oxygen switch co 2 
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Appendix A 

The non-dimensional standard model is 



-AiA Se 



-Gs 



S dA F(S,h,§,t 



dS 

dt ~ 1 - A 1 + B(s + S) ' 1 - A dt ' I- A 



A 7 Se 



-Gs 



(sflpsii(l-5)+flp S ii(s- l))L(A) 



ds 

dt ~ 1 + 5(5 + S) 

+ K 3 ((p - P R )H Ferment (r2 2 -co)- L(A)(sH ¥Sn (l - s) 

s dA 



+ H VS a(s-l))H C Mn(s-S l )) 



A dt 



dp 
dt 



K 2 [L(A)(sH PS1I (l - s) + H PS „(s - l))H Ca i vin (s - Si) 



(p - -PR)H Ferment (^ 2 - CO) 



p dA 
A~dT 



dco 
~di 



= A[K 6 (sH mi (l -s) + H mi (s - 1))L(A) 

- K 5 (coH Rcsv (i2 1 -co) + ^iH Re sp((W - £2i))] 

- V L (co- l)H Loss (co- 1) 



dA 

~dt 



A(R c H G2 (p - yi) + R G P G (p - Yo)Hgi(Y2 - P) 
H G i(p-Yi) -RbH gi (yi -p)) 



dh 
dt 



— ^^Sensivity(^2 _ fw)Hc a l v in(Sl — s) 

[E L S H sL(A) + (1 - E L )P H (p - P R )]L(A) 



where F(S,h,dh/dt,t) is the non-dimensional input func- 
tion and Table II presents standard values for the non- 
dimensional parameters. 
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